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Abstract 

We formulate the algebraic Bethe ansatz solution of the SU{N) vertex models with 
rather general non-diagonal toroidal boundary conditions. The reference states needed 
in the Bethe ansatz construction are found by performing gauge transformations on the 
Boltzmann weights in the manner of Baxter [1]. The structure of the transfer matrix 
eigenvectors consists of multi-particle states over such pseudovacuums and the corre- 
sponding eigenvalues depend crucially on the boundary matrix eigenvalues. We also 
discuss for N = 2 the peculiar case of twisted boundaries associated to singular matrices. 
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1 Introduction 



The study of vertex models have led to important developments in the field of exactly solvable 
models in two dimensions [1]. Their transfer matrices are in general constructed from local 
Boltzmann weights >C^j(A) where A denotes a spectral parameter. This operator can be viewed 
as a matrix on the space of states A representing, for instance, the horizontal degrees of freedom 

L 

of the vertex model on the square lattice. Its matrix elements are operators on ®Vi where 

i=l 

Vi represents the space of states of the vertical degrees of freedom at each site i of a chain of 
size L. The corresponding transfer matrix can be expressed in terms of an ordered product of 
J^AiW operator over the space A denominated monodromy operator X^{\) [2], 

7a(A) = jCal{X)Cal-i{X) . ..Cai{X) (1) 

In terms of the monodromy matrix, a sufficient condition for integrability is the Yang-Baxter 
algebra [2, 3] which reads 

R{x, /x)r^(A) r^(//) = r^(/x) r_^{x)R{x, /x) (2) 

where R{X, /i) is an invertible matrix over complex numbers acting on the tensor product A® A 
space. 

The Yang-Baxter algebra is invariant by the transformation 7^ (A) — > QaXa{X) provided 
that the group of c-numbers matrices satisfies the following property [4] 

[R{X,fj,),gA^QA]^0 (3) 

An immediate consequence of this symmetry is the possibility to define the operator 

T(A) = TrAiGATAiX)] (4) 

which gives origin to generalized families of commuting transfer matrices. 

When the matrix is non singular a quantum spin chain can be associated with the 
transfer matrix (4). For sake of simplicity, consider the usual situation in which the spaces 
A and Vi are isomorphic and that the >C^j(A) is proportional to the exchange operator at 
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certain special point say A = 0. The corresponding one-dimensional Hamiltonian is obtained 
as a logarithmic derivative of the transfer matrix at point A = 0, which reads [4, 11] 

= 2^ -nm — U=o +yL — — |a=o Gl [o) 

i=l 

Clearly, the admissible matrices play the role of more general toroidal boundary condi- 
tions than the particular periodic case when is the identity matrix, the simplest possibility 
satisfying relation (3). Prom the point of view of a vertex model, such general twisted bound- 
ary conditions correspond to the introduction of a line of defects along the infinite direction on 
the cylinder. Though boundary conditions are not expected to influence the inflnite volume 
properties it can change the finite-size corrections of massless systems in a strip of width L 
which contains fundamental informations concerning the underlying conformal field theories 
[5] . For instance, in statistical mechanics boundary conditions provide the means to relate the 
critical behaviour of a variety of different lattice systems such as the Heisenberg spin chain, the 
Ashkin-Teller and the Potts models [6] . In this sense, it is highly desirable to study integrable 
models with as much general boundary conditions as possible. 

If the boundary matrix is diagonal the corresponding transfer matrix (4) can be diago- 
nalized with very little difference from the periodic case because it does not change in a drastic 
way the properties of the monodromy matrix elements. The same does not occur when is 
non-diagonal, starting from the fact that the reference state of the periodic case, essential to 
implement Bethe ansatz approaches, is a priori no longer of utility due to the breaking of the 
original bulk symmetry by the boundary terms. In fact, progress towards solving commuting 
transfer matrices with general twists by Bethe ansatz techniques are modest as compared with 
the literature known for the periodic case, specially for solvable vertex models based on Lie 
algebras, e.g. refs.[7, 8, 9]. To our knowledge, the six vertex model and its higher spin descen- 
dants [10, 11] are the only solvable vertex systems analyzed so far with non-diagonal boundary 
conditions. Even in these cases, the functional relation method used in refs.[10, 11] gives the 
transfer matrix eigenvalues but not information on the corresponding eigenvectors. The latter 
is certainly an important step in the program of solving integrable systems. 
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The purpose of this paper is the formulation of the quantum inverse scattering method for 
the simplest multistate generalization of the six vertex model having N independent degrees 
of freedom on each lattice bond. This turns out to be the isotropic SU (N) vertex model whose 
origin goes back to the work by Uimin [12] and Sutherland [13] on generalized integrable 
Heisenberg chains with higher symmetry. Its corresponding C^iW operators can be written 
as [14, 15] 



where as usual J^j is the identity matrix on the space A^Vi. 

The interesting feature of this system is that the admissible symmetries constitute of arbi- 
trary N X N matrices due to the standard property Pi2^i ® B2 — B2® ^i-Pi2- Therefore 
this provides us a rich variety of possible diagonal and non-diagonal boundary conditions. In 
next section, we present the details of the solution of the eigenvalue problem for the transfer 
matrix in the simplest N = 2 case. Interesting enough, wc find that the Bethe ansatz solu- 
tion depends on the eigenvalue problem related to the boundary matrix. In section 3 we 
generalize these results for arbitrary values of N by using the nested Bethe ansatz approach. 
Our conclusions are presented in section 4 as well as a discussion of singular boundaries for 
the model N = 2. In Appendix A we discuss briefly the completeness of the Hilbert space for 
= 2 and finite L. 

2 Algebraic Bethe ansatz for Heisenberg model 

The purpose of this section is to determine the eigenvalues and the eigenvectors of the following 
transfer matrix 



The operator jC^i{X) is the elementary Boltzmann weights of the isotropic six vertex model 
which can be written as 



(6) 



(7) 



CAi{\) 
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where o"^ j and cr^ j are Pauli matrices acting on the vertical space of states and the weights 
are a(A) = A + 1 and 6(A) = A. The boundary matrix is an arbitrary 2x2 matrix over the 
complex numbers whose matrix elements are denoted by 

\ 



(9) 



5-11 512 
921 922 J 

An essential ingredient of the quantum inverse scattering is the existence of a reference state 
such that the action of the monodromy operator in this state gives as a result a triangular 



matrix. Though each of the operators £^i(A) when acting on the trivial spin up 







or spin down 



states becomes triangular, such property is not extended to the total 



monodromy because the off diagonal elements of Qj^ are in general non-null. Therefore the 
standard ferromagnetic pseudovacuum is not useful when both gi2 and 921 are different from 
zero. In order to find an appropriate reference state we have to introduce a set of gauge 
transformations similar to that used by Baxter [1] in the solution of the eight vertex model. 
We replace the local operators C^i{\) by new matrices >C^j(A) such that [2] 



£^,(A) = Mrl,Cj,^{X)M^ 



(10) 



where Mj are arbitrary invertible 2x2 c-number matrices acting on the space A. After 
performing this gauge transformations the transfer matrix (7) becomes 



r(A) = Trj,[M^^gj,ML+iTj,{\)] 



(11) 



where Tj^{\) = >C^l(A) . . .£^i(A). 

The next step is to look for gauge transformations Mj such that C^j (A) is annihilated for 
instance by its lower left element for arbitrary values of the spectral parameter. Representing 
the matrices Mj by 



M, 



•^3 O 



(12) 



we can conclude [1] that such annihilation property occurs when the ratio is a constant for 



j — 1, . . . , L + 1. As a consequence of that we can choose the local reference state \0)j as 

following that the action of the operator >C_4j(A) in this state is given by 

jC^JX) |o), = ^ I ^ (14) 

where the symbol ^ represents general non-null values. 

The remaining freedom that we have on the matrix elements of Mj is now used to choose 
matrices Mi and M^+i in such way that they transform the boundary matrix into a diagonal 
matrix. More precisely, by imposing that 

M{'g^ML+i = ° (15) 

[O 92 ) 

it follows that the constrains for the first column elements are 

guXL+i + guVL+i = gixi 

g2iXL+i + g22yL+i = giyi (16) 

while for the second column elements we have 

guTL+i + guSL+i = g2ri 

g2irL+l + g22SL+l = g2Sl (17) 

At this point we emphasize our assumption that we are dealing with a non-singular bound- 
ary matrix. While we have an enormous freedom to choose the second column elements the 
same does not occur for the first ones because the ratio ^ needs to be kept fixed to preserve 
triangularity of £^j(A). This latter fact together with relation (16) impose a restriction to this 
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ratio which is precisely the same satisfied by ratio of the components of the eigenvectors of the 
boundary matrix Q^. Therefore, we have two possibihties for the ratio p^^^ = — which are 



{911 - 922) ± J i9u - 922T + 4=912921 . . 

' = -T, (18) 

^921 

Putting now all these informations together it is possible to build up two appropriate global 
reference states |0) by the tensor product 

L I 

io)W = n8) 



p(±) 



(19) 



V 



At this point the state (19) preserves at least the desirable triangular property of the 
total monodromy M^^Qj^ML^xTj^iX). Below we shall see that they are indeed eigenstates of 
the transfer matrix (7) independent of further choices of the elements of the gauge matrices 
My Further progress is made by recasting the Yang-Baxter algebra for the gauge transformed 
monodromy T4(A) in the form of commutation relations for the creation and annihilation fields. 
In order to do that it is convenient to represent Ta(A) by the following 2x2 matrix 



X4(A) 



(20) 



V CiX) D(X) J 

As a consequence of the triangular property (14) we are able to derive important relations 
for the diagonal elements of the transformed monodromy matrix 

A(A)|0)(^) = [a(A)]^^ |0)(^) 

Z5(A)|0)« = lKA)]'»iil-^*^il |0)W pi) 



besides the annihilation property 



C(A) |0)(^^ = (22) 



Now taking into account that gauge matrices Mj are themselves symmetries allowed by the 
property (3) it is not difficult to show that the gauge transformed monodromy Tj^iX) matrix 
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satisfies the same Yang-Baxter algebra as the original monodromy matrix T^{X). In other 
words, we have that Xa{X) satisfies the relation 



(23) 



where in our case the i?-matrix is given by 

( a(A-/i) 

i?(A,//) = 











1 6(A - yu) 
b{X -fx) 1 



(24) 



y U a{\-n) j 

This means that we have the same set of commutation rules of the periodic six vertex model 
[2, 3] however now for the gauged matrix elements. Out of sixteen possible relations three of 
them are of great use, namely 



A{X)B{i^) 
D(A)5(/i) 



a{ii — A) 

b{l^ - A) 

a(A — /i) 



B{^)D{\) - 



b{pi - A) 
1 



B{X)A{,,) 
B{X)D{^i) 



= 



(25) 

(26) 
(27) 



The fields i^(A) are then interpreted as a kind of creation operators over the pseudovacuum 
10)'-''^'* and a natural ansatz for the eigenvectors l^)^^"* of the transfer matrix T(X) is 



10) 



(±) 



n± 



n^(Af)io)(^^ 



(28) 



The eigenvalue problem r(A) 10)^^^ = A(^)(A) 10)^^^ now becomes 

g,A{X)+g2D{X)] \cj>)^^^ = A^^\X) \cj>)^^^ 



(29) 



and it can be solved in the same way as the periodic six vertex model [2], i.e by taking the 
fields A{X) and -D(A) through the creation operators -B(A) with the help of the commutation 
rules (25-27). One peculiarity here, however, is that the calculations involving the action of 
the diagonal fields A{X) and -D(A) over the reference state |0)^^-* requires extra simphfications 



to eliminate unnecessary dependence of the gauge matrices elements. They are carried out by 
using the help of Eqs.(16 - 18) and our final result for the eigenvalues A(='=)(A) are 

AW(A)^ffW[a(A)pn l) , l +9'^'m]'U , 1) ? (30) 

■ - A,; — A — 7J i=l A — A,; + 7J 



1=1 
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provided that the rapidities A^-^^ satisfy the following Bethe ansatz equations 



x(±) 1 



L 
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n /a) .L, . (31) 



9^^^ fA Af ) - Af ) 



where we have performed the convenient shift A^''^'* X^"^^ — ^. The phase factors gf^"*"^ are just 
the eigenvalues of the matrix 



r:. (fi'll + 922) ± V (^11 - 922)^ + 4^12521 

9^ ' = 2 (^^) 

Rather remarkably, we see that the final form of the eigenvalues and Bethe ansatz equations 
resemble much that of the isotropic six vertex model with diagonal boundary if we replace the 
diagonal twists by the eigenvalues of the non-diagonal boundary matrix. The eigenvectors 
can also be thought as multi-particle states in which the integers n± < L play the role of 
particle number sectors. We emphasize, however, that the corresponding basic creation fields 
are much more sophisticated operators than that of the periodic six vertex model [2]. It is 
tempting to think that the two possible ways we have at our disposal to build up the Hilbert 
space is related to the remaining Z2 symmetry allowed by boundary terms. One expects 
therefore that it should be possible to obtain the eigenvalues of the transfer matrix either 
from the |0)*^''''' or |0)*^~^ pseudo vacuums. Indeed, we have verified this fact by numerically 
solving the equations for some values of L and comparing them to exact diagonalization of 
the transfer matrix (7). We note, however, that a given eigenvalue of the transfer matrix is in 
general obtained at different particle sectors n± over the |0)*^^^ reference states. For example, 
the eigenvalue g^[a{\)]^ + 9~[b{X)]^ can be obtained either from the zero-particle state |0)*''''' 
or as a L-particle state over the pseudovacuum |0)^~\ In Appendix A, we present details of our 
study for L = 2 in which Eqs. (30-31) can be solved by analytical means. Our numerical results 
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up to L = 4 suggest that two possible branches of the Bethe ansatz solutions (31) produce the 
complete spectrum of the transfer matrix (7). It would be interesting to further investigate the 
completeness of the Bethe ansatz (31) by adapting the recent arguments developed by Baxter 
[16] to the case of non-diagonal twists. 

We now can derive similar results for the spin chain that commutes with the transfer matrix 
(7). The corresponding spin- 1/2 XXX Hamiltonian follows from expression (5) and it is given 

by 

n 



+1 + C^i CT-i-H + 



with the following boundary condition 



'L+l 



fi'llfi'22 — 912921 



9li 



-gli 



(33) 



-911921 

912922 

\ -'^9ll9l2 ^921922 9ll922 + 912921 J 



\ / + \ 

CTl 



-912 



922 



(34) 



V -1 / 



Its eigenvalues E^^^ = dLogiAi^Hx)] 



U=o are 



n± 



JL 



(35) 



3=1 - 4 ^ 

where satisfy the Bethe ansatz equations (31). 

Our final comment concerns with the comparison between our results (30-31) and that 
of refs.[10, 11] in the isotropic limit case when the trigonometric weights become rational 



functions. We see that they are in accordance for the common non-diagonal boundary Qj^ — 
apart from the fact that our numbers of roots n± can vary up to L while that of 
refs.[10, 11] are fixed at L. This imphes that for non-diagonal boundary conditions the complete 



solution of the isotropic limit does not follows directly from that found for the anisotropic six 
vertex model [10, 11]. This means that even for this particular non-diagonal boundary the 
results (30-31) are novel in the literature. 
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3 Nested Bethe ansatz for SU{N) model 

The purpose of this section is to generahze the results of the previous section for general N. 
We wish to diagonalize the transfer matrix (7), where now the operator jC^i{X) is 

N N N 



a,/3=l a, 13=1 



'a 



(36) 



a=l 



where e-p are the standard Weyl matrices whose elements acting on the space V are [e\j '] 



\kl 



SikSji. In this basis the boundary matrix is generally represented by 



N 

E 



(37) 



As before we have to seek for suitable references states by imposing the gauge transforma- 
tion (10) for each operator (36) and require that they are up triangular when acting on such 



N 



pseudovacuum. Denoting the gauge matrices by Mj = E ''T^ji'^^ P)^ai3 ^^^^ ^^ch 

a,/3=l 

triangular property is fully achieved when the following ratios relations are satisfied 



, a,/3 = l,...,iV 



(38) 



for each j = 1, . . . , L-|- 1. In terms of these ratios the local reference state |0)^. assume the form 



|0), 



P2,l 
PN-1,1 

1 



(39) 



Other important ingredient is the action of the gauge transformed operator £^j(A) over 
the local state of reference. This now can be represented by the following N x N matrix on 
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the space A 



( 



^A, |0), 



V 



a(A)^|0) 








# # ••• # 

■ KA)#r|0), j 



(40) 



ATxAT 



where the variables are given by 

P = 

J a 



N-l 



n 



1 



1,...,N-1 

det[Mj], a^N 



(41) 



Similarly to the previous section we can take advantage of the remaining freedom of the 
elements of the gauge matrices to transform M^^Q^Ml^i into a diagonal matrix. By imposing 
this condition the matrices elements of Mi and M^+i become related by the expression 



N 



7=1 



(42) 



where ga represent the diagonal elements of the transformed boundary matrix. 

Equations (38) and (42) together impose constrains to the possible values ratios Pa,i3 which 
turns out to be same conditions satisfied by the ratio of the components of the eigenvectors of 
boundary matrix Q^. This means that we have N possible choices for / = 1, . . . , A^" and 
consequently from Eq.(39) N kind of suitable local references states |0)^-'^ A natural ansatz 



for the N possible choices of global reference states are 

lo)W = I \®\oy:> I 



n®io) 

3=1 



(43) 



The next step is to write a suitable representation for the gauge transformed monodromy 
matrix in the auxiliary space A. The triangular property (40) suggests us to seek for the 
structure used in nested Bethe ansatz diagonalization of the periodic SU (N) vertex models 
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[14, 15] which is 



r(A) = 



1(A) 
Ci(A) 



^i(A) 



Bn-i{\) 
^iiv-i(A) 



(44) 



\ Cjv-i(A) L>jv-ii(A) ••• Dn-in-iW J 

The triangularity property (40) imphes that the fields Bi{X) play the role of creations 
operators, C'j(A) are annihilation fields while the diagonal operator ^(A) and Dn{\) acts on 
the reference state |0)^'^ as 

1(A) |0)« = [a(A)]^^|0)« (45) 
Ai(A)|0)« 



[KA)] 



/i+l 



l,...,Ar 



(46) 



To construct other eigenvectors we shall use the commutation relations between the gauge 
transformed monodromy matrix elements. The arguments used in section 2 allows to conclude 
that these commutation rules are the same as that already known for the periodic SU{N) 
models [14, 15]. The most useful relations for subsequent derivations are 

a(/i - A) ~ , , T., . 1 



Dij{X)BM 



1 



-B,{^i)A{X) 



b(l^ - A) 
S^(//)A.(A)r«(A-/^)f - 



6(A-/x)"^^"^~"'^'"^' ^'^^ b{X-ii) 

B,{X)Bj{pi) = BMB,{Xy'\X - 



^ Bj{X)DM 



(47) 
(48) 
(49) 



where r^'^\Xyj^ arc the elements of the i?-matrix associated to the SU{N — 1) vertex model. 

In terms of the gauge transformed fields, the eigenvalue problem for the transfer matrix 
T(A) becomes 

r _ N-i "I 

hAiX) + Y: gi+iDu{X) \<l>f^ = A^'\X) \<t>f^ (50) 
1=1 

where l^)*^^"* denotes the eigenvectors. Previous experience with these models [14, 15] suggests 
us to suppose that eigenvectors can be written in terms of the following linear combination 

,(0 _ 5 5 M(M)^-r•"-^•••'^l in\(0 



|0)^'^ = i^a.(Ar)...^.^,(A^f)^M--|o)^ 



(51) 
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where sum over repeated indices a„ = 1, . . . — 1 is assumed. At this stage the components 
j:°-m\-°'^ are thought as coefficients of an arbitrary linear combination that are going to be 
determined a posteriori. 

By carrying on the fields AiX) and -Djj(A) over the multi-particle state (51) wc generate 
terms that are proportional to and those that are not the so-called unwanted terms. The 
first ones will contribute directly to the eigenvalue A('^(A) and are obtained by keeping only 
the first term of the commutation rules (47-48) . These calculations are by now standard in the 
literature and here we present only the main results of the action of the transfer matrix on the 
eigenvector which is 

r(A) 



/i j=i 6(a;- 



A) 



A) 



10) 



(0 



1 



+[KA)]^ n 77- 

i=i KA 
+unwanted terms 



s,,(A^'))...s, ,(A^f)r«(A,{Af''^}),::,7^' 



(52) 



All the pieces entering the above expression can be summarized as follows. The terms 



ai...a I 



T'^^)(A, {A^^'^^D^ 5""^ are transfer matrix elements of an auxiliary inhomogeneous problem 
related to the SU{N — 1) vertex model with twisted boundaries Q defined by 

/ 1 I \ d I ..a I 



T^'\K {>i''}C^1 = r(^HA-A^'))- r(^)(A-A?''))S . . . r^^^A- A^/))r;;^;"^^„.^, (53) 

'"'1 



where Qab denotes the elements of the boundary matrix Q given by 



J3 



\ 















J N 



(54) 



/ N-lxN-l 

The unwanted terms are originated when the variables A^-^'''' of the multi-particle state (51) 
are exchanged with the spectral parameter A due to the second part of the commutation rules 
(47-48). It is possible to collect these terms in closed forms, thanks to the commutation rule 
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(49) which makes possible to relate different ordered multi-particle states. It turns out that 
all the unwanted terms are canceled out provided that the rapidities A^^''^ satisfy the following 
restriction, 



9i 



a(A 



.ai 



(A = Af { Af , ^ = 1, . . . , ^ 



(55) 



Now we reached a point which is fundamental to diagonahze T(-^^(A, {A[ ''}) in order to 
compute the eigenvalues of T(A) and at the same time to solve Eq.(55). This becomes possible 
uire 

AW(A,{Af '')}), namely 



if we require that "'i ^ is an eigenvector of the auxiliary transfer matrix with eigenvalue 



T^'\\ {A^'^Da..^ ^ -i-^ ^ = A^'\X, {Af )})^'^'»i-'^^ (56) 
Inspection of Eq.(52) and Eq.(55) together with Eq.(56) shows that the eigenvalue of T(A) 



IS 



JL, 



I a(AS^'') - A) 



KA)]^n 



1 KAf ''^ - A) 



and the nested Bethe ansatz equations (55) become 



9i 



L+l 



ft 



6(A 



n^A 



(1,0 _ ^(1,0 



a(A 



(1,0 _ ^(1,0) 



6(A 



(1,0 _ ^(1,0^ 



=1 b{X - A 



AW(A = Ar'^{A(^'')}), i^l,...,m[ 

(58) 



In order to solve the eigenvalue problem (56) it is necessary to introduce a second algebraic 
Bethe ansatz for the eigenvectors ^""'i" "^ Because the boundary matrix § is triangular 
there is no need to perform gauge transformations to find an appropriate reference state for 
T^^\X, {Aj-^''^}). We can use, for instance, the usual ferromagnetic pseudovacuum build up by 



tensor product of elementary {N — 1) -dimensional 



/1\ 




vectors. As a result the solution 



VO/iv-i 
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(56) becomes very similar to that of the periodic SU{N — 1) vertex model in the presence 
of inhomogeneities. Since this problem has been extensively discussed in the literature we 
will only present our final results for the main eigenvalue problem (50). It turns out that the 
eigenvalues of the transfer matrix A^''*(A) is given by 



Al'l(A;{AS'''>},...,{Ar''"}) = SiTfeWA^n^ 



h j=i b{X) - A) 



^-2 /I < a(A - Af '')) "Ji^ a(Af - A) 



+[KA)]"E^^+i^n : L) n 



.^/^^7.YiS="i - )) U\ 6(Af+^'')-A) 



+ii>(m^i&T n (59) 

Jn j=i b{X-X) ") 

The rapidities {A-''''''} k = 1, . . . ,iV parameterize the multi-particle states of the nesting 
problem at step k and are required to satisfy the following nested Bethe ansatz equations 



91 fur' 



^ a{)^^ - Af '^) < ai>^^ - Af '')) 

\}. a{>^^ - Af '^)) M 5(Af - Af'^)) ^'"^ 



/^/^ "^.^ a(Af - Af-^'^') < a(Af - Af V^U(Af - Af ''') 



(61) 



KA^^-Af-^'O a(Af - Af-'^) 6(Af+^''^-A 



A; = 2, . . . , AT - 2 



^iv /^^l/^ y 6(aS^-^''^ - Af -^'')) " y a{\f-'''^ - Af-^-')) ^ ^ 

The final step is to carry out simplifications on the phase factors Qiju+i with the help of 
the constraints (38) and (42). After a cumbersome algebra it is possible to show that such 
factors are just the eigenvalues ^f*^*^ of the boundary matrix Qj\^. To make sure that the different 
possibilities we have at hand for the ratios p^^^ do not lead us to singular gauge matrices we 
choose to order them for each Z-th choice of pseudovacuum by ^{0^^(/+iV)for; = l,...,Ar. 
Taking into account this ordering as well as performing the shifts {A^'' '''} — > {A^''''^} — k/2 our 
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result (59) for the eigenvalue becomes 

j=i Aj- - A - 2 



Ar-2 K \ _ \(fc,0 I fc+2 i,(fc+l,0 X I 1-fc 



(iV-1,0 , AT+l 



A- Ar^ + 



^=1 A-Ar"'''^ + 



(Ar-1,0 _L Af-i 
'i 



and the nested Bethe ansatz equations can be compactly written as 



\{a,l) 6a,i 
Aj + — 



A 



(«.0 <5o,i 



L 



9" ' TT TT \ ^fc + ~ . _ 1 I, 



nn /w) ^M) cly i-h...,m^ a = l,...,7V-l (64) 



2~J ^ b=l k=i \ - Ajt 2 



where Caf, is the Cartan matrix elements of the SU (N) Lie algebra. 

Wc sec that the results for the eigenvalues and the Bethe ansatz equations is similar to 
that expected from the SU{N) vertex model with diagonal twists giving by the eigenvalues 
of the boundary matrix Q^. It remains to be investigated whether this interesting feature is 
particular of the SU{N) symmetry or it also works in other isotropic vertex models such as 
those invariant by the 0{N) and Sp{2N) Lie algebras. 



4 Concluding remarks 

In this paper we have been able to apply the quantum inverse scattering program to solve 
exactly the isotropic SU (iV) vertex model with non-diagonal twisted boundary conditions. We 
find that the eigenvectors can be constructed in terms of multi-particle states over N possible 
pseudovacuums. The Bethe ansatz results for the eigenvalues are similar to that of the SU {N) 
model with diagonal boundaries in which the eigenvalues of the boundary matrix play the 
role of the diagonal twists. 

We expect that our results can be generalized without further difficulties to accommodate 
the solution of vertex models based on the SL{N\M) super Lie algebra [13, 17, 18] with 
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general non-diagonal twists. These will include interesting systems of correlated electrons 
on a lattice such as the one-dimensional supersymmetric t-J model [19] and the so-called 
Essler, Korepin and Schoutens superconducting model [20] with arbitrary symmetry breaking 
boundary conditions. With more effort we hope that our approach can be further generalized 
to include the trigonometric deformation of those vertex models based on the Ug[SL{N\M)] 
symmetry. In these cases, however, we recall that the possible matrices compatible with 
integrability belong to a smaller group formed by one-dimensional dilatations and the discrete 
Zn+m symmetry. 

Other interesting issue that deserves investigation is the situation when the boundary ma- 
trix becomes singular. For example, one would like to ask it is still possible to exhibit 
eigenvectors of the transfer matrix that are given by direct tensor product of A^-dimensional 
vectors such as the reference states of sections 2 and 3. We have studied this problem in the 
simplest case N — 2 and surprisingly we found a family of such states which are 

n /_a22\ L /^\ 

i0)(") = n® n ® h n=o,i,...,L (65) 

whose corresponding eigenvalues A'^")(A) have also the following simple factorized form 

A^-\\) = {gu + g22M\t--[b{\)r (66) 

This result prompted us to study further properties of the transfer matrix (7) when 
is a singular matrix. Our study for finite L up to six sites reveals that the roots of the 
characteristic polynomial of T(A) are exactly the eigenvalues (66) whose degeneracy is the 
binomial coefficient (i„ — {L-n)\n\ • singular boundary matrix T(A) becomes 

defective since it has fewer than 2^ independent eigenvectors. To each distinct eigenvalue 
A(")(A) we find only one eigenvector which is precisely the state (65) and therefore the total 
number of independent states is L -|- 1. These results are strong evidences that T{X) behaves 
as a non-derogatory matrix and we conjecture that its Jordan decomposition for arbitrary L 
should be 

T{X)^diag{Jo,Ji,---,JL) (67) 
17 



where Jn is a d„ x dn Jordan matrix is given by 

( A(")(A) 1 

Jn 



A('^)(A) 1 








(68) 



V ••• A^-^A); 

This turns out to be a remarkable example how boundary conditions can change in a drastic 
way the Hilbert space of integrable models. At this point it is natural to ask what happens 
to the Bethe ansatz states (28) when one gradually varies the boundary matrix towards the 
singular manifold. In particular, if we can figure out the kind of Bethe states in each sector 
n± that should collapse to the eigenvectors (65). A precise answer to this question as well as 
possible generalizations of these results for arbitrary N has eluded us so far. 



Acknowledgements 

We would like to thank the organizers of the Workshop Flux, Charge, Topology and Statistics 
2003, Amsterdam, where part of this work was carried out. The authors G.A.P Ribeiro and W. 
Galleas thank Fapesp (Fundagao de Amparo a Pesquisa do Estado de Sao Paulo) for financial 
support. The work of M.J. Martins has been partially support by the Brazilian Research 
Council-CNPq and Fapesp. 



Appendix A : Completeness for L = 2 

This appendix is concerned with the study of the completeness of the Bethe ansatz solution (31) 
for L — 2, i.e. that all four eigenvalues of the transfer matrix are obtained either by starting 
with lO)*-^-* or with |0)^~''. Let us first begin with lO)*-^-* whose corresponding eigenvalue Ao^''(A) 
is clearly 

Al^\X) ^ g^^\a{X)r + g(-^[b{X)r (A.l) 
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The next step is to solve the Bethe ansatz equations for the one-particle state -B(Ai) |0)*^''"^ 
As a result we find two possible rapidities given by 



giving us the following one-particle Af eigenvalues 

A± = a{\)h{\) + g^-^) ± V^W^^^ (A.3) 
Repeating similar exercise for the two-particle state ^(Ai)^(A2) |0)*^^'* we have 

j(+) + ;,I-)±2/v/9'+V-' 
- 2fa.-.-,W) 

and the corresponding eigenvalue is 

k,{X) ^ g^-\a{X)f + g^^\b{X)f (A.5) 

An exact diagonahzation of the transfer matrix (7) corroborates these four possible eigen- 
values for L = 2. Note also that (A.5) is exactly the eigenvalue associated to reference state 
|0)''~\ The others eigenvalues (A.l) and (A.3) are easily obtained from |0)~ by noticing that 
to each solution A-"^^ one can find the corresponding A-~^ through the refiection Aj-~^ = — A-"*"^ 
symmetry. We have also investigated numerically this problem for L — 3, 4 and found that 
both references states can lead to the complete spectrum of T(A). 
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